function[pa,qa,pb,qb]=pdebc(xa,ua,xb,ub,t)

pa = [ua(1);ua(2)];qa=[0;0];

pb = [ub(1);ub(2)];qb=[0;0];
